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Abstract 

Purpose Geospatial details about land use are necessary to 
assess its potential impacts on biodiversity. Geographic 
information systems (GIS) are adept at modeling land use 
in a spatially explicit manner, while life cycle assessment 
(LCA) does not conventionally utilize geospatial information. 
This study presents a proof-of-concept approach for coupling 
GIS and LCA for biodiversity assessments of land use and 
applies it to a case study of ethanol production from agricultural 
crops in California. 

Materials and methods In Part 2 of this paper series, 
four biodiversity impact indicators are presented and 
discussed, which use the inventory data on habitat 
composition and sizes from the GIS-based inventory 
modeling in Part 1 (Geyer et al. 2010). The concepts used 
to develop characterization models are hemeroby, species 
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richness, species abundance, and species evenness. The 
biodiversity assessments based on species richness, abun¬ 
dance, and evenness use a species-habitat suitability matrix 
which relates 443 terrestrial vertebrate species native to 
California to the 29 habitat types that occur in the study area. 
Results and discussion The structural similarities and 
differences of all four characterization models are discussed 
in some detail. Characterization factors and indicator results 
are calculated for each of the four characterization models 
and the 11 different land use scenarios from Part 1 of this 
paper series. For the sugar beet production scenarios, the 
indicator results are in fairly good agreement. For the com 
production scenarios, however, they come to fundamentally 
different results. The overall approach of using GIS-based 
inventory data on land use together with information on 
species-habitat relationships is not only feasible but also 
grounded in ecological science and well connected with 
existing life cycle impact assessment efforts. 

Conclusions Excluding biodiversity impacts from land use 
significantly limits the scope of LCA. Accounting for land 
use in inventory modeling is dramatically enhanced if LCA 
is coupled with GIS. The resulting inventory data are a 
sound basis for biodiversity impact assessments, in partic¬ 
ular if coupled with information on species-habitat relation¬ 
ships. However, much more case studies and structural 
analysis of indicators is required, together with an evalu¬ 
ation framework that enables comparisons and ranking of 
indicators. 

Keywords Abundance • Biodiversity impacts • Bioethanol • 
Evenness • Geographic information systems (GIS) • 
GIS-based inventory modeling • Habitats • Hemeroby • 
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1 Introduction and background 

The purpose of life cycle assessment is to synthesize a 
comprehensive environmental profile of products (Rebitzer 
et al. 2004). For many products, the area of land required to 
produce the raw materials, to manufacture a functional unit 
of the product, and to dispose of it at the end of life is 
negligible and is appropriately ignored in life cycle 
assessment (LCA). Some production processes, however, 
transform and occupy substantial areas of land and thereby 
generate significant environmental impact. Agricultural 
production, including for renewable fuels, is a particularly 
noteworthy example, especially if it transforms high 
biodiversity lands such as tropical forest into lower 
diversity cropland (Searchinger et al. 2008). Concern that 
ignoring land use will yield incomplete environmental 
profiles has spawned a growing literature on ideas for 
measuring land use impacts within the LCA methodology 
(Baitz 1998; Lindeijer 2000; Schenck 2001; Brentrup et al. 
2002; Mila i Canals et al. 2006; Koellner and Scholz 2007; 
Michelsen 2008). A UNEP/SETAC committee has recently 
proposed a conceptual framework for land use in LCA 
(Mila i Canals et al. 2007). Biodiversity is explicitly 
included as one aspect of land quality in this framework. 
How exactly to incorporate the framework into standard 
LCA practice is still a very active area of research, 
however. 

There are several serious obstacles to achieving this 
goal. First, land use is an inherently spatial activity. The 
geographic variation in the productivity of soils means that 
the land area required to produce a liter of bioethanol 
depends on where the feedstock is grown. In addition, the 
impact on biodiversity of transforming land depends not 
only on the new use but also on the habitat type and its 
complement of species that would be lost through that 
transformation. Therefore, the impact is a function of place, 
not just of units of product. LCA, however, is conventionally 
aspatial, mostly relying on average values. Overcoming this 
obstacle will require more creativity than merely adding 
another impact category linked to the reference flow of the 
product or a simple tally of land area transformed and 
occupied. 

This relationship of impact and place implies a second 
obstacle to land use impact assessment. The cumulative 
impact on biodiversity will not scale linearly with volume 
of production. Each marginal unit of production will 
transform and occupy a varying marginal land area 
(because of variation in productivity), and each additional 
unit of land may be different in its initial habitat and 
diversity than another unit or the “average.” LCA method¬ 
ology is typically linear, where elementary flows and 
characterization models scale linearly with activity levels 
and economic output. 


Overcoming these substantial technical obstacles still 
leaves a fundamental issue of how to measure as complex a 
concept as biodiversity. Simply put, biodiversity is the 
diversity of life as a whole. The Convention on Biological 
Diversity defined it as “diversity within species, between 
species and of ecosystems” (UNEP 1992). Diversity 
between species is the most common way that biodiversity 
is interpreted in practice, such as in biological conservation 
planning (Margules and Pressey 2000). But which species? 
Any location is likely to contain thousands of species. Only 
the distributions of vertebrates (e.g., birds, mammals) are 
reasonably well known for large regions. Plant distributions 
are also known for many parts of the world, but typically at 
a coarse spatial resolution, such as a district or county or in 
Europe by arbitrary 50x50-km squares (Jalas and Suominen 
1972-1994). Consequently, decisions must be made about 
which species will be practical to include in an impact 
assessment. Another practical decision regards which aspects 
of diversity to measure. The most basic metric is the number 
of species or species richness (Magurran 2004). Trans¬ 
forming 1 ha of land may make the site unsuitable for a 
given species, reducing the richness of that site or its alpha 
diversity (Whittaker 1972). What is of greater concern for 
conservation is maintaining the total richness of an eco- 
region, or gamma diversity, which is a product of the alpha 
diversity of its different communities and the dissimilarity 
(beta diversity) between them (Whittaker 1972). To diminish 
gamma diversity would require the transformation of all 
suitable habitats for a species in an ecoregion to eliminate it 
from the regional species pool. Therefore, gamma diversity 
or richness may not be a very sensitive impact indicator. 
Ecologists often incorporate measures of species abundance 
(absolute size of their populations) or evenness (relative size 
of populations) along with richness for more nuanced 
measures of biodiversity. As a result, many different 
measures of biodiversity have been proposed (Magurran 
2004). Koellner and Scholz (2008) offered an LCA 
biodiversity impact indicator based on species richness. 
Scholes and Biggs (2005) proposed a “biodiversity intactness 
indicator” that estimates the average relative abundance of all 
plant and vertebrate species compared to a reference condi¬ 
tion. Burke et al. (2008) developed a form of abundance- 
based indicator called the “biotope method” that accounted 
for the relative abundance (or rarity) of narrowly distributed 
plants. 

Geyer et al. (2010) in Part 1 presented a life cycle 
inventory method that used a geographic information 
system to calculate elementary flows of habitat types. The 
objective of this article was to demonstrate how a suite of 
characterization models can be derived analytically from 
two basic pieces of information: the elementary flow data 
and knowledge of habitat preferences of wildlife species 
encapsulated in a habitat suitability matrix. Four biodiversity 
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impact indicators, and their associated vectors of character¬ 
ization factors for habitat types, are presented. The indicators 
range from a simple method based on overall degradation 
to habitat types, to a richness-based index, to richness 
plus either abundance or evenness. Potential biodiversity 
impacts were then assessed for three bioethanol feedstocks at 
different production levels in the southern part of California’s 
San Joaquin Valley, USA (Fig. 1 in Geyer et al. 2010). This 
region is renowned both for its agriculture (Umbach 2002) 
and imperiled biodiversity (USFWS 1998). We do not 
conclude that any of these indicators should become the 
standard for life cycle impact assessment. Rather, we seek to 
demonstrate that a wealth of information on biodiversity can 
be mined from relatively basic spatial and ecological 
data with LCA-compatible methodology. We also explore 
whether this LCA-compatible methodology can accommo¬ 
date nonlinear biodiversity impacts in response to production 
level. 

2 Materials and methods 

2.1 The scope of the impact assessment 

The first part of this paper series explained how GIS-based 
inventory modeling was employed to generate ten plausible 
fuel crop production scenarios in the four southernmost 
counties of California’s San Joaquin Valley (Geyer et al. 
2010). The scenarios included three potential bioethanol 
feedstocks—com grain, com grain and stover, and sugar 
beets—at various production levels. The production levels 
were set in order to meet increasing percentages of 
projected demand for ethanol in California in 2010. Each 
scenario consists of a detailed land use map of the entire 
study region, with the required amount of com and sugar 
beet fields displacing current land use. A map of current 
land use is added as a baseline scenario. To be able to 
assess biodiversity impacts of the fuel crop production 
scenarios, the inventory model calculated habitat composi¬ 
tion vectors a from the land use maps by identifying each 
land use type as a habitat type from a list of 59 mutually 
exclusive and collectively exhaustive standard habitat types 
for California (Mayer and Laudenslayer 1988). Twenty- 
nine of the 59 different habitat types are present in the study 
area. Com fields, for example, belong to the irrigated grain 
crops (IGR) habitat type, while sugar beet fields belong to 
irrigated row and field crops (IRF). The areas of all land use 
types belonging to the same habitat type are summed up. 
The baseline scenario and each fuel crop production 
scenario is thus characterized by a vector of elementary 
input flows a,j , j e [1-29], which contains the areas (in 
hectares) of 29 habitat types that together make up the total 

study area, A = ^ aj = 2,538,482 ha. 

j 


The second part of this paper series explores ways in 
which these elementary flows can be used to assess the 
impacts that these fuel crop production scenarios may have on 
biodiversity. The remainder of this paper is thus dedicated to 
deriving and discussing characterization models and factors 
that translate the habitat composition vectors a into four 
different biodiversity indicators. They are derived from the 
following four concepts: 

• Hemeroby, or naturalness, as proposed by Brentmp et 
al. (2002) 

• Species richness, as proposed by Koellner and Scholz 
(2008) 

• Species richness and abundance, developed by the 
authors 

• Species richness and evenness, based on the Simpson 
index (Simpson 1949) 

Hemeroby is a property of the land use or habitat type 
itself, which means that an indicator based on hemeroby 
does not require any information about the species the 
habitat supports. The other three indicators, on the 
other hand, need to relate the land or habitat type to 
the species that are found in it, which is the major challenge 
of biodiversity impact assessment models. To achieve 
this, three of the four indicators introduced in this 
paper rely heavily on a species-habitat suitability matrix, 
HS G R 443 x R 29 , which characterizes the suitability of 
each of the 29 habitats for 443 terrestrial vertebrate species 
native to the study area in California. Each matrix element 
is either 0 or 1, with hsy= 1 meaning that habitat j is suitable 
for species /, i.e., its requirements for food, cover, and 
reproduction are met, and hsy =0 that it is not. The species- 
habitat suitability matrix HS is based on data from the 
California Wildlife Habitat Relationships (CWHR) system 
(CDFG 2005). The CWHR information system used expert 
opinion to rate the suitability of each of the 59 habitat types 
for each native vertebrate species in four ordinal classes— 
nil, low, moderate, and high. For this study, we treated the 
moderate and high ratings as suitable habitats, i.e., 
converted them to hsy= 1 for species i in habitat y, and the 
nil and low ratings as unsuitable habitats, i.e., converted 
them to hsy= 0 for species i in habitat y. It should be noted 
that specific habitat patches are not necessarily inhabited by 
all the species; they are suitable for and therefore should be 
understood as potential habitats. 

2.2 Impact assessment based on the hemeroby concept 

To assess the potential impact of a given habitat composi¬ 
tion a on biodiversity, models are required that relate land 
use to habitats and habitats to biodiversity. One very basic 
model is the hemeroby concept which, according to 
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Kowarik (1999), measures “human influence on ecosys¬ 
tems” by determining its “deviation from naturalness as a 
result of specific land use types” (Brentrup et al. 2002). 
According to Brentrup et al. (2002), the concept has been 
introduced in ecology 55 years ago, and its use in LCA has 
been proposed by Klopffer and Renner 1995. Brentrup et 
al. (2002) converted the original cardinal scale of increasing 
human influence into 11 evenly spaced characterization 
values that increase from 0 to 1 in 0.1 increments. The 
hemeroby concept has been specifically proposed as an 
alternative to other methods that assess land use impacts on 
species diversity, such as the one proposed by Koellner 
(2003). Hemeroby-based impact assessment does not 
require information on species diversity, such as the 
species-habitat suitability matrix HS used in the other three 
characterization models discussed in this paper. Instead, it is 
based on the hypothesis that land use impacts on 
biodiversity increase monotonically with the level of human 
influence and interference. For this study, we adopted the 
characterization factors proposed by Brentmp et al. (2002) for 
European land cover/use types, called naturalness degrada¬ 
tion potential, and assigned them to the corresponding 
habitat types in our study area (CF1 in Table 1). 

The resulting impact assessment follows the structure of 
classification and characterization as described in ISO 
14044 (ISO 2006). First, the elements cij of the habitat 
composition vector a are identified as the elementary flows 
that contribute to the biodiversity impact category. Next, 
the hemeroby-based biodiversity impact potential BIP1 is 
calculated by multiplying the elementary flows aj with their 
characterization factors CFlye [0,1]: 

BIPl = NCFl/ x a r (1) 

j 

It can be seen that BIPle[0, A\, with B IP 1=^4 for 
maximum potential impact on biodiversity and BIP1=0 for 
minimum potential impact on biodiversity, e.g., if the entire 
study area were lacustrine or riverine habitat. While a 
hemeroby-based biodiversity impact assessment is very 
easy to implement, it has little theoretical or empirical 
underpinning. While a relationship between human influ¬ 
ence and biodiversity is plausible, the hemeroby concept 
offers no environmental mechanism between the LCI 
results (habitat composition vector a) and the category 
endpoint, which is lack of biodiversity. Such environmental 
mechanism can be based on theoretical models, empirical 
data, or combinations of the two. The CWHR system, for 
example, is a comprehensive data resource containing life 
history, geographic range, habitat relationships, and man¬ 
agement information on 694 species of amphibians, 
reptiles, birds, and mammals known to occur in the state 
(CDFG 2005). It has been used to derive the species- 


habitat suitability matrix HS introduced in the previous 
section. The remaining three characterization models 
illustrate different ways in which such a matrix can be 
used for biodiversity impact assessment. 

2.3 Impact assessment based on species richness 

The most basic way to use the species-habitat suitability 
matrix HS for biodiversity impact assessment is to simply 
count the number of species for which a habitat type is 
deemed suitable, i.e., to sum up the matrix elements of each 
column. Each habitat type j is then characterized by the 
number of species Sj, for which it is potential habitat: 

Sj = '£ hs V ( 2 ) 

i 

The use of species richness data has been suggested in 
Koellner and Scholz (2008) who proposed a set of 
characterization factors called ecosystem damage potentials. 
Following their linear approach to convert species richness 
data into characterization factors (CF ; ) for biodiversity impact 
assessment yields: 

CF2 ; = 1 - A = 1 - -L ^ hSir (3) 

Oref O re f j 

The value for the reference species richness, S rQ f, used in 
this paper is 443, which is the number of native terrestrial 
vertebrates for which at least one of the occurring 29 habitat 
types is suitable. It can be seen from the above equation 
that CF2 / e[0, 1]. The value of the characterization factor 
decreases as the number of species increases for which the 
habitat type j is suitable. This follows the convention that 
characterization factors measure the potential adverse 
impact that an elementary flow has for the underlying 
environmental concern. The resulting biodiversity impact 
potential BIP2 for a given land use scenario is now 
calculated as the scalar product of the vector of character¬ 
ization factors and the habitat composition vector: 

BIP2 = CF2y x cij = A — —— h s ijdj with 

j SvQf j i 

S Ye f = 443 and^4 = aj = 2,538,482 ha. 

j 

(4) 

It can be seen from the equation above that BIP2e [0, A], 
with BIP2=^4 if the habitat composition of the study area 
provides no potential habitat to any species and BIP2=0 if 
it provides potential habitat to all species everywhere. 
Realistic values will be somewhere in between these 
extremes. 
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Table 1 Habitat types and their characterization factors 





Habitat type (je[ 1-29]) 

Hemeroby 

(CFly) 

No. of suited 
species (S 7 ) 

Pure richness 
(CF2y) 

Richness and 
abundance (CF3 7 ) 

Richness and evenness 
(CF4, x io 9 ) a 

Agricultural and urban habitats 

DOR—Deciduous orchard 

0.7 

99 

0.78 

0.90 

1.60 

EOR—Evergreen orchard 

0.7 

43 

0.90 

0.97 

0.93 

IGR—Irrigated grain crops 

0.8 

103 

0.77 

0.92 

2.03 

IRF—Irrigated row and field crops 

0.8 

65 

0.85 

0.97 

1.56 

IRH—Irrigated hayfield 

0.8 

164 

0.63 

0.83 

2.70 

PAS—Pasture 

0.55 

57 

0.87 

0.94 

0.97 

RIC—Rice 

0.8 

96 

0.78 

0.89 

1.47 

URB—Urban 

0.95 

153 

0.65 

0.80 

1.96 

YIN—Vineyard 

0.7 

62 

0.86 

0.96 

1.30 

Native habitats 

AGS—Annual grassland 

0.1 

163 

0.63 

0.84 

2.77 

ASC—Alkali desert scmb 

0.1 

78 

0.82 

0.89 

0.84 

BAR—Barren 

0.8 

85 

0.81 

0.87 

0.80 

BOP—Blue oak-foothill pine 

0.1 

206 

0.53 

0.71 

2.34 

BOW—Blue oak woodland 

0.1 

205 

0.54 

0.71 

2.37 

COW—Coastal oak woodland 

0.1 

201 

0.55 

0.72 

2.33 

CRC—Chamise-redshank chaparral 

0.1 

115 

0.74 

0.84 

1.39 

CSC—Coastal scmb 

0.1 

145 

0.67 

0.81 

1.82 

DSC—Desert scmb 

0.1 

96 

0.78 

0.84 

0.81 

EU C—Eucalyptus 

0.7 

86 

0.81 

0.90 

1.30 

FEW—Freshwater emergent wetland 

0.1 

136 

0.69 

0.82 

1.65 

JUN—Juniper 

0.1 

118 

0.73 

0.84 

1.35 

LAC—Lacustrine 

0 

120 

0.73 

0.81 

1.10 

MCH—Mixed chaparral 

0.1 

154 

0.65 

0.77 

1.64 

MHW—Montane hardwood 

0.1 

178 

0.60 

0.73 

1.74 

PGS—Perennial grassland 

0.1 

160 

0.64 

0.84 

2.64 

RIV—Riverine 

0 

110 

0.75 

0.83 

1.04 

SGB—Sagebmsh 

0.1 

123 

0.72 

0.83 

1.48 

VOW—Valley oak woodland 

0.1 

203 

0.54 

0.72 

2.33 

VRI—Valley-foothill riparian 

0.1 

230 

0.48 

0.66 

2.35 


Com is part of the IGR habitat type; sugar beets are part of IRF. CF4y are shown here for the baseline habitat a b only 
a With this characterization model, every scenario has its own set of characterization factors CF4 7 


Koellner and Scholz (2008) also suggest a nonlinear 
approach which uses the natural logarithm of the normal¬ 
ized species richness ln(*S)/A re f) rather than normalized 
species richness itself to derive characterization factors. 
The idea behind this is to have a characterization model that 
reflects the redundant species hypothesis, which says that 
the marginal utility of ecosystem services is diminishing 
with increasing species richness. Since the presented 
research is concerned with biodiversity rather than ecosys¬ 
tem service assessment, the linear approach has been used. 

While species richness is probably the most prominent 
aspect of biodiversity, it is only one part of a very complex 
concept that encapsulates the full range of life, including 
diversity within species, between species, and of ecosys¬ 


tems (UNEP 1992). Other important aspects of diversity 
between species, also called alpha diversity, are species 
abundance and evenness, i.e., the number of individuals for 
each present species and the similarity of the population sizes 
of similar species at the same trophic level. The character¬ 
isation models shown above and proposed by Koellner and 
Scholz (2008) yield pure species richness indicators which 
reflect neither species abundance nor evenness. 

2.4 Impact assessment based on species richness 
and habitat frequency 

In order to go beyond pure species richness, characterization 
models for biodiversity impact assessment need to capture 
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additional aspects of biodiversity, such as abundance and 
evenness. One natural extension of the approach discussed in 
the previous section is to weight each species according to its 
rarity, i.e., lack of abundance (Lindner 2008). Using r t as a 
yet-to-be-defmed rarity weight, the characterization factors 
can now be written as: 


CF3/ = 1-V Vi x hsu. 

J s K i v 



Rarity rather than abundance is used in the equation 
above since characterization factors conventionally increase 
with increasing environmental impact. The rarer the species 
supported by habitat type y, the lower should thus be the 
habitat’s characterization factor. Data on species abundance 
or rarity is usually generated through sampling and is even 
more demanding in terms of sampling techniques and 
effort. Sampling is not very practical in the context of LCA 
and even impossible if future land use changes are forecast. 
However, the species-habitat suitability matrix HS can be 
used to derive information related to relative species 
abundance, assuming population size increases with area. 
Multiplying a given habitat composition vector a with HS 
yields a vector that contains the total area of suitable habitat 
hi (in hectares) for each of the covered species, which are 
443 native terrestrial vertebrates in the presented research: 


h = HS x a with h t = hsy x aj , i G [1 — 443]. (6) 

j 


abundance or rarity of species is a property of a given 
region at a given point in time. It is therefore not fixed but 
rather a function of time and also the chosen study area. 
Habitat frequency therefore has to be defined for a given 
region at a given time. In this paper, the baseline land use 
scenario, i.e., the study area in its current habitat compo- 
sition a, has been used to derive habitat frequency data h t . 
In principle, rarity weights could be also obtained other¬ 
wise, e.g., by using threatened species lists or information 
on species’ global habitat and population sizes. 

Given the baseline habitat composition a b , r t =0 if the 
entire study area A is potential habitat for species /, and r t = 
1 if species i has no potential habitat in the study area A. 
Using the above definition for species rarity weights r h the 
resulting characterization factors are calculated as: 



In this paper, A ref =443 and ^4=2,538,482 ha. It can be 
seen that the species-habitat suitability matrix HS is now 
used twice: once to derive potential species richness Sj for 
each habitat type a 7 - and again to derive species rarity 
weights r t from species’ habitat frequencies hi. The 
characterization factors can also be written as: 


It can be seen that for each species z, the above equation 
sums up the areas of all suitable habitat types, i.e., those aj 
for which hsy= 1. Given a habitat composition a, the 
resulting h t of a species z can be used as a proxy for its 
abundance. To obtain unitless values between zero and one, 

/V 

habitat frequencies h t are calculated by expressing each 
total area of suitable habitat h t as a proportion of the total 
study area: 

hi = ^ with A = aj and hi G [0,1]. (7) 

J 


CF3 ' = 1 - E' hs «+ ^ (- hSij x a J) xhSij - 

( 10 ) 

In this form, it can be seen that these characterization 
factors are the sum of the pure species richness factors and 
a correction term that accounts for habitat frequency and 
thus reflects species abundance. Using matrix notation to 
rewrite the correction term for habitat frequency, the 
resulting biodiversity impact potential BIP3 for a given 
land use scenario is: 


Since habitats support many different species, J] h t » A 

_ /\ 

and therefore 22 hi >> 1. Rarity weights r t can now be 
defined in a number of ways, the simplest of which is: 

n = 1 - hi with r; G [0,1]. (8) 


j 

+ 



1 

BIP3 = > CF3/ x aj = A - 

^ S re f 

rei J i 

(HS x a b ) r HS x a. 


hsy x aj 


1 


5ref X A 



One consequence of using the total area of potential 
habitat h t in order to account for habitat frequency is the 
fact that characterization factors that are based on these 
parameters become functions of a given habitat composi¬ 
tion vector a. This is simply a reflection of the fact that 


It should be noted that there are various other functions 

that could be used to derive rarity weights from the habitat 

/\ 

frequency /z z , including a simple reciprocal transformation 

/\ 

Vi = l/hi as well as various forms of logistic functions, 

/V 

such as — l/{a + exp(/z z - — b)). At a future point, it 
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would be useful to compare these different options in terms 

of their usefulness for biodiversity impact assessment. In 

/\ 

this paper, the linear transformation r t = 1 — hi has been 
chosen since it is the most straightforward and transparent 
option and thus an appropriate starting point. 

2.5 Impact assessment based on species richness 
and habitat evenness 

The third and last way in which we employed the species- 
habitat suitability matrix HS for biodiversity impact assess¬ 
ment is based on the Simpson index, which is a widely used 
biodiversity indicator that accounts for species richness and 
evenness (Simpson 1949). For a given sample of N 
individuals, the Simpson index is calculated as: 

SI = N(N- 1) ^2 n ‘ x ~ 9’ with N = ZN (12) 

with rii being the number of individuals in the sample that 
belong to species /. Simpson (1949) showed that the index 
calculates the best estimate, given the sample size, for the 
probability that two individuals randomly picked from the 
study area belong to the same species. The indicator is 
sensitive to species richness and evenness since its value 
decreases with increasing number of species and increasing 
evenness of species distribution. Lower index values 
therefore indicate higher biodiversity. Since sampling is not 
very practical in the context of LCA, the species-habitat 
suitability matrix HS is employed to derive a biodiversity 
indicator that has the same mathematical structure as the 
Simpson index but is based on total potential habitat per 
species h t rather than number of sampled individuals «/. This 
indicator is calculated as follows: 

BIP4 = (hf) 2 with hi — hsy x cij and H — h t . 

i j i 

(13) 

Like the Simpson index, the indicator value decreases 
with increasing number of species and increasing evenness 
of potential habitat distribution. Noting that X] (hi) 2 = h T h 
and h = HS x d, the indicator can be written in matrix 
notation as: 


BIP4 



(HS x <J) r HS x a. 



Comparing Eqs. 11 and 14, it can be seen that BIP4 
bears strong resemblance to the habitat frequency correc¬ 
tion term of indicator BIP3 of the previous section. This is 


not entirely surprising since habitat frequency and evenness 
are related concepts. The two differences are that the 
Simpson index-based indicator has a different normaliza¬ 
tion factor and uses the habitat composition of the assessed 
land use scenario a rather than that of the baseline scenario 
a b to assess species abundance. The characterization factors 
of the Simpson index-based indicator can be derived from 

BIP = ^2 CFj x a,j and are 

j 

CF4 / - = — —y (hsy x cij) xhsy. 

(E E hs a x aA ' 1 

(15) 

It can be seen that the elementary flows dj are required to 
calculate their characterization factors, which is not the case 
in conventional life cycle impact assessment where charac¬ 
terization factors tend to be independent of the elementary 
flows. However, ISO 14044 only states that the conversion 
of LCI results into indicator results uses characterization 
factors and does not require that characterization factors 
be constants. Like the characterization factors from 
Sections 2.2 and 2.3, CF1 and CF2, the Simpson index- 
based characterization factors CF4 do not require a baseline 
land use scenario. Unlike the characterization factors CF1 
and CF2, but like the characterization factor CF3 from 
Section 2.4, the Simpson index-based characterization 
factors CF4 are a function of habitat composition vectors. A 
comprehensive theoretical analysis of the structural similari¬ 
ties and differences between all four discussed indicators was 
outside the scope of the presented research and thus left to 
future research. Instead, the four characterization models were 
applied to the inventory modeling results from the fuel crop 
production scenarios developed and discussed in Geyer et al. 
(2010). The remainder of this paper discusses the results of 
these biodiversity impact assessments. 


3 Results 

The characterization factors of all four characterization 
models are shown in Table 1. It can be seen that the 
hemeroby-based factors CF1 offer the least amount of 
differentiation between the different habitat types. Of the 
nine agricultural or urban habitats, seven have a value of 
either 0.7 or 0.8, while 18 of the 20 native habitats have a 
value of 0 or 0.1. Note that the CF1 for com and sugar 
beets are both 0.8. The only difference in impacts of the 
two crops for a given volume of ethanol will therefore be 
due to the difference in their yields and the habitats they 
replace. CF2 simply reflects the number of species for 
which each habitat is deemed unsuitable. As would be 
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expected, the average value for the agricultural or urban 
habitats (0.79) is higher than the average value for the 
native habitats (0.67). Comparing CF1 and CF2 shows that 
species richness is not always well reflected through the 
hemeroby concept. Urban habitats and irrigated hayfields 
are just as species-rich as mixed chaparral and annual 
grassland even though their hemeroby values are dramat¬ 
ically different. Weighting each species based on its rarity 
increases the characterization factors for every habitat since 
CF3y > CF2y <^> J2 r ih s y <J2hsy and avg(/y) = 0.74. The 
rarity weights have a more dramatic impact on habitats 
with large Sj = ^ hsy, which narrows the range of the CF3 y 
to 0.32 from 0.42 for the CF2. Interestingly, value range 
and standard deviation for the group of agricultural or urban 
habitats decrease, as do value range and standard deviation 
for the group of native habitats. This means that accounting 
for species rarity makes the agricultural or urban habitats 
more similar to each other, which is also true for the native 
habitats. The Simpson index-based characterization factors 
CF4 cannot be easily compared with the other three factors. 
Whereas species richness and abundance are calculated for 
each individual species in isolation, evenness relates all 443 
species to each other. As a result, characterization factors 
CF4 are different for each land use scenario, which makes 
their interpretation difficult. It is more meaningful to 
discuss BIP4 directly rather than CF4. 

Table 2 contains the indicator results for all 11 land use 
scenarios and four characterization models. The first three 
indicators, BIP1, BIP2 and BIP3, are essentially weighted 
sums of the habitat areas aj describing a land use scenario. 


Table 2 Biodiversity impact potentials for each scenario and 
characterization model 


Scenario 

BIP1 

BIP2 

BIP3 

BIP4 

Baseline 

1,422,653 

1,883,105 

2,265,559 

0.005102 

CR10 

1,429,939 

1,880,993 

2,265,702 

0.005134 

CR25 

1,441,063 

1,876,994 

2,265,355 

0.005193 

CS10 

1,428,327 

1,881,487 

2,265,698 

0.005127 

CS25 

1,436,638 

1,878,588 

2,265,510 

0.005168 

CS50 

1,462,724 

1,876,109 

2,266,862 

0.005288 

SB 10 

1,425,056 

1,884,985 

2,267,159 

0.005125 

SB25 

1,428,488 

1,887,760 

2,269,494 

0.005159 

SB50 

1,434,963 

1,891,901 

2,273,072 

0.005218 

SB75 

1,508,484 

1,917,307 

2,289,524 

0.005431 

SB 100 

1,707,205 

1,984,195 

2,333,725 

0.006327 


The goal has been defined as 1.38 billion kilograms of ethanol. The 
numbers in the scenario names indicate the percentage of goal 
achievement (for example, CS25 means 25% of ethanol production 
goal achieved with com including stover conversion) 

CR com grains only, CS com grains and stover, SB sugar beets 


Since the weights (characterization factors) are all between 
0 and 1, the resulting indicator values are all between 0 and 
A, the total study area. Smaller indicator values indicate 
smaller adverse biodiversity impacts. Comparing the 
absolute indicator results across the first three characteriza¬ 
tion models, however, is not very meaningful. To facilitate 
comparisons across all four indicators, all indicator results 
have been normalized in the following way: 


, . _ BIP(X) - BIP (baseline) 

v y BIP(^100) - BIP (baseline) 



This means that BIP (baseline) = 0 and BIP (SB 100) = 1. 
The results can be seen in Table 3 and have also been 
visualized separately for each ethanol feedstock in Figs. 1, 2, 
and 3. 

It is important to discuss the results in the context of the 
underlying inventory modeling, which is discussed in detail 
in Part 1 of this paper series (Geyer et al. 2010). While the 
entire study region with its 2,538,482 ha and 29 habitat 
types has been modeled for each fuel crop production 
scenario, the vast majority of habitat type area changes 
happens within just six habitat types, which are deciduous 
orchard (DOR), IGR, IRF, irrigated hayfield (IRH), 
vineyard (VIN), and annual grassland (AGS). Table 4 
shows the percentage to which the fuel crop habitat 
increases come from the habitat types listed above. 
A«dor/^^igr(CR 10) = —28%, for example, means that 
28% of the land converted to com fields in the CR10 
scenario was of the habitat type DOR. It can be seen that 
com production scenarios convert overwhelmingly DOR, 
IRF, IRH, and VIN into IGR. The first three sugar beet 


Table 3 Percentage change of biodiversity impact potentials relative 
to baseline scenario for each scenario and characterization model 


Scenario 

BIP1 (%) 

BIP2 (%) 

BIP3 (%) 

BIP4 (%) 

Baseline 

0.0 

0.0 

0.0 

0.0 

CR10 

2.6 

-2.1 

0.2 

2.6 

CR25 

6.5 

-6.0 

-0.3 

7.4 

CS10 

2.0 

-1.6 

0.2 

2.0 

CS25 

4.9 

-4.5 

-0.1 

5.4 

CS50 

14.1 

-6.9 

1.9 

15.2 

SB 10 

0.8 

1.9 

2.3 

1.8 

SB25 

2.1 

4.6 

5.8 

4.6 

SB50 

4.3 

8.7 

11.0 

9.5 

SB75 

30.2 

33.8 

35.2 

26.8 

SB 100 

100.0 

100.0 

100.0 

100.0 


Percentage changes have been normalized to scenario SB 100, which 
has the highest biodiversity impact potential for each characterization 
model 


4^ Springer 













700 


Int J Life Cycle Assess (2010) 15:692-703 



Fig. 1 Percentage change of biodiversity impact potentials relative to 
baseline scenario for com grain only production scenarios ( CR ) 


Fig. 3 Percentage change of biodiversity impact potentials relative to 
baseline scenario for sugar beet production scenarios (SB) 


production scenarios, SB 10, SB25, and SB50, convert 
overwhelmingly DOR, IGR, IRH, and VIN into IRF. In 
the remaining two sugar beet production scenarios, SB75 
and SB 100, large areas of AGS are turned into sugar beet 
fields as well. The relationship between the characterization 
factors of these six habitat types is thus the key to 
understanding the indicator results. 

Figures 1 and 2 show the biodiversity indicator results 
for the com production scenarios (with and without stover). 
While the hemeroby- and Simpson index-based indicators 
BIP1 and BIP4 are in good agreement, the other two 
indicators come to very different conclusions. According to 
BIP1 and BIP4, increasing corn production increases 
potential impact on biodiversity. In contrast, BIP2 indicates 
a reduction in biodiversity impacts as com production 
increases, while biodiversity impacts according to BIP3 are 
roughly constant. In other words, the indicators fundamen¬ 
tally disagree on the biodiversity impacts of the com 
production scenarios. The discrepancy in the indicator 
results follows from the relative differences in the charac¬ 
terization factors. CFl/Gi? > CFl^A: E [DOR, IRF, IRH, 
VIN], so converting a mix of DOR, IRF, IRH, and VIN 
into IGR will increase BIP1. In contrast, C¥2 1GR < CVl^k 
E [DOR, IRF, VIN], which explains why BIP2 decreases 



Fig. 2 Percentage change of biodiversity impact potentials relative to 
baseline scenario for com grain and stover production scenarios ( CS) 


with increasing com production. CF3 IGR is larger than the 
values for DOR and IRH, but smaller than those for IRF 
and VIN, which explains why BIP3 is hardly affected by 
the com production scenarios. CF4 jgr^) is considerably 
larger than the values for DOR, IRF, and VIN, which 
together make up 70-80% of the converted habitats. 
However, the Simpson index-based characterization factors 
change with each crop production scenario, which compli¬ 
cates this kind of analysis. 

Figure 3 shows the biodiversity indicator results for the 
sugar beet production scenarios. Here, the results for all 
four indicators are in fairly good agreement. All four 
indicators conclude that increasing sugar beet production in 
the study area increases the potential impact on biodiver¬ 
sity. For all four indicators, the increase in potential 
biodiversity impact is modest for the first three production 
targets, SB 10, SB25, and SB50, but rather dramatic for 
SB75 and SB 100. The most significant reason for this 
nonlinearity is that SB75 and SB 100 convert large areas of 
AGS into IRF (see Table 4). However, they also increas¬ 
ingly use other native habitat types such as ASC, BOW, and 
MCH (see Geyer et al. 2010 for habitat composition details) 
with characterization factors much lower than those for IRF. 

Tables 2 and 3 also allow comparisons across fuel crop 
types for given ethanol output and indicator choice. It can 
be seen that the four indicators come to conflicting results, 
which is due to their different assessments of the com 
production scenarios. BIP1 and BIP4 conclude that for a 
given ethanol target, sugar beet generates less potential 
biodiversity impact than com, with or without stover. BIP2 
and BIP3 come to the opposite conclusion. 


4 Discussion 

This paper discusses and applies four different character¬ 
ization models for an important environmental mechanism, 
which is the potential impact of land use on biodiversity. 
It is important to recognize that other environmental 
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Table 4 Habitat area changes m 
percent of fuel crop habitat 

Scenario 

DOR (%) 

IGR (%) 

IRF (%) 

IRH (%) 

VIN (%) 

AGS (%) 

Other (%) 

increase, i.e., Ao/A^gr for 









com production scenarios and 

CR10 

-28 


-25 

-16 

-24 

-3 

-4 

A< 2 /A< 2 IRF for sugar beet 

CR25 

-24 


-24 

-16 

-30 

-3 

-3 

production scenarios 

CS10 

-29 


-25 

-16 

-23 

-3 

-4 


CS25 

-24 


-25 

-17 

-28 

-2 

-4 


CS50 

-22 


-25 

-16 

-24 

-6 

-7 


SB 10 

-26 

-20 


-11 

-37 

-1 

-5 


SB25 

-23 

-19 


-12 

-41 

-1 

-4 


SB50 

-21 

-21 


-10 

-43 

-1 

-4 


SB75 

-16 

-15 


-6 

-27 

-30 

-6 


SB 100 

-9 

-10 


-4 

-14 

-49 

-14 


interventions, such as toxic releases, can also impact 
biodiversity, while land use can impact safeguard objectives 
other than biodiversity, such as water purification and 
nutrient cycling. Such mechanisms are outside the scope of 
the presented research. How to generate inventory data that 
describe land use of product systems and contain enough 
pertinent information to characterize its potential impact on 
biodiversity was presented in the first part of this paper 
series. This second part focuses entirely on the character¬ 
ization of the inventory results, i.e., the habitat composition 
vectors a. 

Apart from the very simple hemeroby-based approach, 
all discussed characterization models are based on the 
species-habitat suitability matrix HS. This is motivated in 
equal measure by the widespread availability of this type of 
information and its usefulness in characterizing biodiversi¬ 
ty. HS is very similar to a presence-absence matrix which, 
according to Arita et al. (2008), is the basic analytical tool 
underlying many biogeographic species diversity studies 
and concepts. Arita et al. (2008) emphasize the many 
different ways in which a presence-absence matrix can be 
analyzed. As with the HS matrix, analysis by columns 
alone yields species richness per site, whereas analysis by 
row yields species distribution. Considering the interaction 
of rows and columns offers more sophisticated insights into 
the relationship between biogeography and biodiversity. 
The matrix approach presented here and in Arita et al. 
(2008) appears to be an excellent way to link LCA with 
ecology and should be studied in more detail in theory and 
practice. For example, the matrix elements of HS were 
reduced in this study to either 0 or 1. They could also 
reflect the four suitability levels of the CWHR information 
system, e.g. 0, 0.3, 0.7, and 1 for no, low, moderate, and 
high suitability, respectively, or be modeled as continuous 
values (Gottschalk et al. 2005). It is not obvious how this 
would change the indicator results since the matrix values 
would increase for some species-habitat relationships and 
decrease for others. Species-habitat relationships have been 
compiled around the world, particularly in the USA through 


the Gap Analysis Program (Scott et al. 1993), but are still 
relatively scarce in tropical regions (Gottschalk et al. 2005). 

The emerging standard in the assessment of land use 
impacts within LCA is to describe land use in terms of area 
affected, duration of the land use, and change in land 
quality due to the land use (Mila i Canals et al. 2007). It is 
also customary to describe land use as a combination of 
transformation and occupation impacts (Mila i Canals et al. 
2007). Transformation is typically defined as a more or less 
instantaneous change in the properties of the land (mea¬ 
sured in area). Occupation is usually defined as the actual 
use of the land for some purpose (measured in area-time) 
where changes in the properties of the land are accidental 
and tolerated rather than intended. Standard land use 
scenarios consist of a transformation process followed by 
a time-limited occupation and subsequent relaxation. A lot 
of attention has been paid to this temporal aspect of land 
use (see, e.g., Guinee 2002; Koellner and Scholz 2007; 
Mila i Canals et al. 2007). In order to focus entirely on the 
unresolved problem of how to measure “land quality” in 
terms of biodiversity impact, the presented research is 
based on transformations followed by infinite occupation. 
At a later point, the emerging characterization models 
should be extended to include the temporal aspects of how 
land use impacts biodiversity. 

A much discussed issue for characterization models of 
land use impacts is the need for a reference state against 
which the impacts are to be measured (Mila i Canals et al. 
2007). It should be noted that three out of the four 
discussed characterization models do not require a refer¬ 
ence state to calculate a biodiversity indicator result for a 
given habitat composition vector. Hemeroby and potential 
species richness are absolute qualities of each habitat type, 
while the habitat evenness calculated by indicator BIP4 is 
an absolute property of the given habitat composition. Only 
in the case of BIP3 was the calculation of potential species 
abundance h based on the baseline habitat composition 
vector, which therefore constituted a reference state for 
species abundance. In principle, however, the rarity weights 
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Ti — 1 — hi used for BIP3 could be based on the studied 
habitat composition a rather than the baseline a b . Typically, 
the biodiversity impacts of habitat composition vectors will 
be compared against each other. The object of study is then 
a land use change, with the initial land use scenario as a de 
facto reference state. Future research should further look 
into the relationship between reference states and attribu- 
tional versus consequential perspectives. 

For all three crop types, com grains only (CR), com 
grains and stover (CS), and sugar beets (SB), and all four 
indicators, the relationship between ethanol output and 
potential biodiversity impact is distinctly nonlinear. This 
nonlinearity therefore appears to be a robust effect 
independent of fuel crop and indicator choice. There are 
two factors that cause this nonlinearity: The first factor is 
the nonlinear relationship between fuel crop area and 
ethanol output (see Fig. 3 in Geyer et al. 2010). As ethanol 
output increases, less and less productive areas are 
converted for fuel crop production, which causes the 
ethanol area function to be convex. While the convexity is 
a property of the profit-based land conversion modeling, 
the nonlinearity is generic. The second factor is the fact that 
increasing ethanol output leads to conversion of increasing 
amounts of native habitats, which tend to have smaller 
characterization factors than agricultural habitats. This 
increases the potential biodiversity impact per hectare of 
converted land. That both nonlinearities are being 
accounted for is due to the use of total ethanol output as 
reference flow rather than some fixed arbitrary amount. It is 
straightforward to derive average and marginal data from 
total input-output relationships, but not vice versa. This 
suggests that inventory and impact models that are based on 
absolute activity or output levels are better suited to capture 
nonlinearities than models that use average or marginal 
data. 

One conspicuous outcome of the San Joaquin Valley 
case study is that even simple indicators based on different 
aspects of biodiversity can come to fundamentally different 
results. To better understand the reasons for this, more 
practical applications and systematic theoretical analysis of 
the four indicators are required. In retrospect, the presented 
case study can be seen as challenging since it mostly 
consists of conversions from one agricultural land use/ 
habitat type to another. It is possible that more dramatic 
conversions from exclusively natural habitats to exclusively 
agricultural habitats would produce more consistent results. 
However, in the long mn, it would be desirable to develop 
characterization models that are capable of meaningful and 
robust biodiversity assessments of relatively subtle land use 
changes. This includes the ability to distinguish not only 
between different agricultural land uses, such as orchards 
versus grain crops, but also practices, such organic versus 
conventional farming or different tillage practices. 


Finally, the use of terrestrial vertebrates as indicator 
species needs to be reexamined in future research. It is 
currently unclear how robust characterization factors and 
impact indicator results are with regard to the selected species 
group. It is known that species richness among taxonomic 
groups is only weakly correlated (Wolters et al. 2006). This 
issue will be a major aspect of our subsequent research. 

5 Conclusions 

Comprehensive inventory data on habitat types and areas, 
derived from GIS-based inventory modeling of land use, 
are a sound basis from which to conduct biodiversity 
impact assessment of land use. The next step is to develop 
characterization factors that relate such inventory data to 
the category indicator, which is biodiversity. While charac¬ 
terization factors for different land use/habitat types can be 
derived directly from concepts such as hemeroby, more 
detailed environmental mechanisms, such as theoretical or 
empirical species-habitat relationships, provide more detailed 
and scientific characterization models. Species-habitat suit¬ 
ability matrices, for example, are powerful tools that allow the 
derivation of data on potential species richness, abundance, 
and evenness. The use of species-habitat suitability matrices 
for biodiversity impact assessment of land use is very 
promising. The presented case study of fuel crop production 
scenarios in California’s San Joaquin Valley, however, shows 
that at least in the case of habitat changes limited mostly to 
agricultural types, different characterization models generate 
contradicting results, even if they are based on the same 
matrix. 

As can be seen from the equations, the biodiversity 
indicators BIP2, BIP3, and BIP4 use only the species- 
habitat suitability matrix HS and the habitat composition 
vectors from the inventory analysis. The GIS-based 
inventory modeling generates comprehensive land use and 
habitat maps of the entire study area and therefore contains 
considerably more information than the habitat composition 
vectors convey. This additional information could, in 
principle, be used for more advanced biodiversity assess¬ 
ments. When aggregating all habitat patches of the same 
habitat type, for example, patches could be eliminated if 
they are below a defined minimum size and therefore not 
viable for given species. In addition to total area, the 
connectivity or fragmentation of habitats could be assessed. 
Another possibility would be to model and account for edge 
effects. The potential for using spatial data generated by 
GIS-based inventory modeling for biodiversity impact 
assessment, as well as other impact categories, is enormous 
and offers a wealth of research opportunities. 

There are currently no methods available to assess how 
suitable, meaningful, and scientific a given biodiversity 
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indicator is, which is why we refrained from making 
recommendations. There is thus an urgent need for an 
evaluation framework that allows for comparative assess¬ 
ments of the strengths and weaknesses of these and other 
potential biodiversity indicators of land use. Such a 
framework should identify and prioritize the biodiversity 
aspects that indicators should reflect and also specify how 
proposed indicators can be measured against these criteria 
and compared with each other. Otherwise, we will face 
further proliferation of competing characterization models. 
Comparative assessments should be equally based on 
practical applications of the studied characterization models 
and their conceptual and systematic theoretical analysis. 
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